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ABSTRACT 



Context. Theoretical predictions from models of red giant branch stars are a valuable tool for various applications in astrophysics 
ranging from galactic chemical evolution to studies of exoplanetary systems. 

Aims. We use the radiative transfer code Optim3D and realistic 3D radiative-hydrodynamical (RHD) surface convection simulations 
of red giants to explore the impact of granulation on interferometric observables. We assess how 3D simulations of surface convection 
can be validated against observations. 

Methods. We compute intensity maps for the 3D simulation snapshots in two filters: in the optical at 5000 ± 300 A and in the K band 
2.14 ± 0.26 yum FLUOR filter, corresponding to the wavelength-range of instruments mounted on the CHARA interferometer. From 
the intensity maps, we construct images of the stellar disks, accounting for center-to-limb variations. We then derive interferometric 
visibility amplitudes and phases. We study their behavior with position angle and wavelength, and compare them with CHARA 
observations of the red giant star HD 214868. 

Results. We provide average limb-darkening coefficients for different metallicities and wavelength-ranges. We detail the prospects 
for the detection and characterization of granulation and center-to-limb variations of red giant stars with today's interferometers. 
Regarding interferometric observables, we find that the effect of convective-related surface structures depends on metallicity and 
surface gravity. We provided theoretical closure phases that should be incorporated into the analysis of red giant planet companion 
closure phase signals. We estimate 3D- ID corrections to stellar radii determination: 3D models are ~ 3.5% smaller to ~ 1% larger in 
the optical with respect to ID, and roughly 0.5 to 1.5% smaller in the infrared. Even if these corrections are small, they are important 
to properly set the zero point of effective temperature scale derived by interferometry and to strengthen the confidence of existing red 
giant catalogues of calibrating stars for interferometry. Finally, we show that our RHD simulations provide an excellent fit to the red 
giant HD 214868 even though more observations are needed at higher spatial frequencies and shorter wavelength. 
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1. Introduction 



stars: atmospheres - hydrodynamics - radiative transfer - techniques: interferometric stars: 



Red giant branch stars have evolved from the main sequence 
and are powered by hydrogen burning in a thin shell surround- 
ing t heir helium core. T heir masses are typically less than ~2.0 
M G (ISalaris et al. 2002), their effective temperatures range from 
~4 000 to ~5 100 K, depending on metallicity, their surface gray - 
ities from ~3.5 to ~ 1.0 in log g dHek ker & Meiendez l2007h. and 
their r adii from ~3 to -70 R Q (Ivan Belle et al.ll999HBaines et al.l 
l2010l) . The determination of fundamental parameters of red gi- 
ant stars is of great importance in astrophysics: (i) red giant stars 
will be used as tracers of the morphology and evolution of the 
Galaxy in the framework o f the Gaia mission dPerrvman et al.l 
l200lHLindegren et al.|[2008h . as they are intrinsically bright and 
thus can probe regions obscured by interstellar extinction; (ii) 
they are extensively used for spectroscopic ele mental abundance 
analyses of distant ste l lar populations (e.g. iMc William et al.l 
1 19951: iRvan et al.ll 19961; lFulbrighdl2000l); (iii) they are relat ively 
easy targets in open dYong et al.l 20051: I Carnev et al.ll2005l) and 
globular clusters (see iGratton et al 1 12004 for a review) and can 
be used for a number of purposes, including measuring Galactic 



metallicity gradients; (i v) red giant stars have been included in 
planet search surve ys dSato et al.1 12001 iDollinger et al.1 120071; 
I Johnson etai1l2008l) . 

Three-dimensional hydrodynamical modeling of convec- 
tion of stellar surfaces has been already extensively em- 
ployed to study the effects of photospheric inhomogeneities 
and velocity fields on the formation of spectral lines in 
a number of cases, including the Sun, dwarfs and sub- 
giants (e.g. 



. [A splund et al 
l Asplund et al.l 
ISbordone et al 



...1999t I Asplund & Garcia Perezll2001 
20091: ICaffau et al.l 120101: iBehara et al 
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fl) and red giants (e.g. Collet et all 



2010; 



2007 
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(e.g. ICott 

120091; IWende et all 120091: [Kucinskas et alj |2009l) . Such Simula 
tions are paramount for an accurate quantitative analysis of ob- 
served data. These studies provide valuable information to prop- 
erly understand astrophysical processes such as stellar nucle- 
osynthesis, internal mixing mechanisms in stars, and Galactic 
chemical evolution. 

Within this framework, the correct interpretation of interfer- 
ometric observables is very important to properly recover fun- 
damental stellar parameters via a correct description of limb- 
darkening laws. Once angular diameters are known, the com- 
bination with bolometric fluxes and parallaxes provides effec- 
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tive temperatures (T e ff) and physical radii. A proper chracteriza- 
tion of these quantities is also relevant from the point of view of 
stellar evolution, since red giants are populating the coolest and 
most luminous part of the Hertzsprung-Russell (HR) diagram. 

In this work, we present interferometric predictions obtained 
from three-dimensional surface convection simulations of red gi- 
ant stars with different effective temperatures, surface gravities, 
and metallicities. 



2. Surface convection simulations and radiative 
transfer calculations 

2.1. Radiative-hydrodynamical simulations of red giant stars 

We adopt here time-dependent, three-dimensional (3D), 
radiative-hydrodynamical (RHD) sur face convect i on sim ula- 
tions of red giant stars ge nerated by Collet et alJ (I2007I) and 
ICollet et all (l2009t) with the lStein & Nordlundl (Il998h code and 
with the Stagger-CodeO respectively. In these simulations, the 
equations of conservation of mass, momentum, and energy are 
solved for representative rectangular volumes located across the 
stars' optical surfaces. Heat losses and gains due to radiation 
are accounted for by solving the 3D radiative transfer equation 
in the local thermodynamic equilibrium (LTE) approximation. 
The simulations employ realistic input p hysics: the equation o f 
state is an updated version of the one bv lMihalas et al. (jl988l). 
Continuous and line opacitie s are taken from Gustafsso n et al.l 
dl975b ; iKuruczl (1 1 992l fl 9931) . The simulation domains are cho- 
sen large enough to cover at least ten pressure scale heights ver- 
tically and to allow for about ten granules to develop at the sur- 
face; they are periodic horizontally, while their top and bottom 
boundaries are left open. The parameters of the adopted simu- 
lations are given in Table Q] For all simulations, we assum e a 
solar chemical composition from lGrevesse & Sauvalld 19981) but 
with all metal abundances scaled proportionally to the desired 
iron abundance. The simulations assume a constant gravitational 
acceleration parallel to the vertical axis; hence, there is a degen- 
eracy between stellar mass and radius (g oc M/R 2 ). Nonetheless, 
in the present work, we assume for each simulation a mass of 
0.8 M - typical of red giants with stellar parameters like the 
ones considered here - to estimate the spatial frequency scale in 
Section|5] 



2.2. Post processing radiative transfer code 

We used the 3D pur e-LTE radiative transfer code Optim3D 
(IChiavassa et al.l2009l) to compute intensity maps from the snap- 
shots of the RHD simulations listed in Table Q] The code takes 
into account the Doppler shifts due to convective motions. The 
radiative transfer equation is solved monochromatically using 
extinction pre-tabulated coefficients as a function of tempera- 
ture, density, and wavelength. The lookup tables were computed 
for the same chemical compositions as the RHD simulations us- 
ing the same extensive atomic and mol ecular opacity data as th e 
latest generation of MARCS models dGustafsson et al.l l2008t) . 
We assumed a zero micro-turbulence since the velocity fields 
inherent in 3D models are expected to self-consistently and ad- 
equately account for non -thermal Doppler broadening of spec- 
tral lines (lAsplundll2000t) . The temperature and density ranges 
spanned by the tables are optimized for the values encountered 
in the RHD simulations. 

Optim3D computes the emerging intensities for vertical rays cast 
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Fig. 1. Intensity maps in the optical top-hat filter of Fig. [3] for 
the simulation with logg =1.6 and [Fe/H]=-3.0 (TableQ]) with 
H—l.QO, 0.81, 0.58, and 0.40. The intensities typically ranges 
from3xl0 4 to 1.2 x 10 6 ergcrrr 2 s~' A -1 . 



through the computational b ox, for all requi r ed wa velengths, 
with the method described by IChiavassa et al ] (120091) . The pro- 
cedure is repeated after tilting the computational box by an an- 
gle with respect to the line of sight (vertical axis) and rotat- 
ing it azimuthally by an angle <f>. The final result is a spatially 
resolved intensity spectrum at different angles. This implemen- 
tation of Optm3D enables radiative transfer calculations with 
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Table 1. Simulations of red giant stars used in this work. Note that the mass has been fixed to determine the radius. 



3 





[Fe/H] 


logg 


x, y, z-dimensions 


x, y, z-dimensions 


M 


R, 


[K] 




[cgs] 


[Mm] 


[grid points] 


[M e ] 


[Ro] 


4697±18^, 


0.0 


2.2 


1250x1250x610 


100x100x125 


0.8 


12.9 


4717±12~ 


-1.0 


2.2 


1125x1125x415 


100x100x125 


0.8 


12.9 


5035+13 6 


-2.0 


2.2 


1150x1150x430 


100x100x125 


0.8 


12.9 


5128±10 6 


-3.0 


2.2 


1150x1150x430 


100x100x125 


0.8 


12.9 


4627±14 c 


-3.0 


1.6 


3700x3700x1100 


480x480x240 


0.8 


23.0 



" Temporal average and standard deviation of the emergent effective temperatures 
* ICollet et al.l ( 120071) 

c ICollet et al.l d2009h and Collet et al. 2010, in prep. 



RH P simulations eith er in the star-in-a-box configuration (e.g. 
Chiavassa et al 1 l2009h or in the box-in-a-star configuration (this 
work). 

3. Sphere tiling: from local images to a semi-global 
intensity maps 

In this Section, we aim to obtain an image of the stellar disk as a 
nearby observer would see it as required to extract the interfero- 
metric observables. 

The computational domain of each simulation represents 
only a small portion of the stellar surface (see Section [27Tb . To 
overcome this limitation, and at the same time account for limb- 
darkening effects, we computed surface intensity maps for dif- 
ferent inclinations with respect to the vertical (0-angles) and for 
representative series of simulation snapshot^ and used them to 
tile a spherical surface. The actual value of the 0-angle used to 
generate each map depended on the position (longitude and lat- 
itude) of the tile on the sphere and was determined by comput- 
ing the angular distance of the tile from the center of the stel- 
lar disk (0° in longitude and latitude on the sphere). For prac- 
tical reasons, we actually computed intensity maps for a set of 
15 predefined values of //= cos(<9)=[1.000, 0.989, 0.978, 0.946, 
0.913, 0.861, 0.809, 0.739, 0.669, 0.584, 0.500, 0.404, 0.309, 
0.206, 0. 104] (Fig.[TJ and interpolated linearly in between them. 
Mapping onto a spherical surface induces distortions especially 
at high latitudes and longitudes; we accounted for that by appro- 
priately cropping the square-shaped intensity maps when defin- 
ing the spherical tiles. 

In addition, we took into consideration statistical tile-to-tile 
fluctuations in the number of granules and in their shape and size 
and computed intensity maps by selecting snapshots at random 
from each simulation's time-series. This random selection pro- 
cess also avoided that the assumption of periodic boundary con- 
ditions resulted in a tiled spherical surface globally displaying an 
artifactual periodic granulation pattern. Based on the stellar radii 
estimates given in Table [TJ and on the sizes of the simulations' 
domains, we required 14 to 23 tiles to cover half a great circle 
from side to side on the sphere, depending on the surface grav- 
ity of the models (log g- 1.6 and 2.2, respectively). The adopted 
number of /v-angles for the pre-computed intensity maps is there- 
fore sufficiently large to ensure an accurate representation of the 
center-to-limb variations for the present work. To produce the 
final stellar disk images, we performed an orthographic projec- 
tion of the tiled spheres on a plane perpendicular to the line-of- 
sight (6 = 0°). The orthographic projection returned images of 

2 For each simulation, we typically selected ~30 snapshots, taken at 
regular intervals and covering ~100 - 250 hours of stellar time 



the globes in which distortions are greatest toward the rim of 
the hemisphere where distances are compressed. The sketch in 
Fig.|2]outlines our method we used to tile spherical surfaces and 
constructing the synthetic stellar disk images. 
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Fig. 2. Toy model representing our method used to tile spherical 
surfaces and constructing the synthetic stellar disk images. 



For the computation of the interferometric observables, we 
consider two spectral bands: (i) in the optical between 4700 
and 5300 A and (ii) in the infrared at 2.14+0.26 /im FLUOR 
filter, Fig. [3] These wavelength intervals have been chosen be- 
cause they correspond to the range of instruments m ounted 
on CHARA interferometer (Iten Brumm elaar et al. 1 12001 . These 
telescopes are particularly suited for observations of red giant 
stars because of its long baselines allow to resolve objects with 
small apparent diameters. The resulting intensity maps reported 

are normalized to the filter transmission as: ^^Trr- where h is 

J T{A)dA 

the intensity and T (A) is the transmission curve of the filter at 
a certain wavelength. The spectra in Fig.[3]have been computed 
along rays of four yu-angles [0.88, 0.65, 0.55, 0.34] and four </>- 
angles [0°, 90°, 180°, 270°], after which we performed a disk 
integration and a temporal average over all selected snapshots. 
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Wavelength {/j.m) Wavelength (fj.m) 

Fig. 3. Left panel: synthetic spectra for different 3D model red giant atmospheres in the wavelength range 4700-5300 A . The 
intensity maps corresponding to this filter have been com puted assuming a top-hat filter (dashed line). This wavelength range 
corresponds to the VEGA instrument dMourard et al.ll2009l) mounted on CHARA with a spectral resolution A/AA — 5000. From the 
bottom: the black curve refers to the simulation with logg = 2.2 and [Fe/H]=0.0 and no offset applied, the red curve to logg = 2.2 
and [Fe/H]=-1.0 with an offset of 1, the green curve to logg = 2.2 and [Fe/H]=-2.0 with an offset of 2, the light blue curve to 
logg = 2.2 and [Fe/H]=-3.0 with an offset of 3, and the dark blue curve logg = 1.6 and |~Fe/ H"|=-3.0 with an off set of 4 (see 
Table [TJ. Right panel: transmission curve (dashed line) of the FLUOR (2.14 ± 0.26 pm) filter dMerand et al.l 12006b mounted on 
CHARA. Colors and offsets are the same as in left panel. 



4. Three-dimensional limb darkening coefficients 

In this Section we derive average center-to-limb intensity pro- 
files with their limb darkening (LD) laws for the simulations in 
Table Q] Figure |4] shows the resulting synthetic stellar disk im- 
ages for the [Fe/H]=-3.0 and logg=1.6 simulation in both fil- 
ters of Fig. [3] The surface of the simulated giants stars appear 
irregular, permeated with structures. The center-to-limb varia- 
tions are visible in both images, but are more pronounced in the 
optical filter. This is due to the integrated extinction by atomic 
and molecular lines in the upper layers which is more impor- 
tant in the optical range than in the FLUOR filter: this difference 
further increases for metal-poor models (see Fig. [3}. 

We then derived azimuthally averaged intensity profiles for 
every synthetic stellar disk image from the simulations, with two 
examples shown in Fig. [5] The profiles were constructed using 
rings regularly spaced in p — cos(8) for p < 1 (i.e. r/R* < 1), 
with 6 the angle between the line of sight and the radial direc- 
tion. The parameter p is related to the impact parameter r/R* 
through r/R* = yT- p 2 , where R* is the stellar radius reported 
in Table Q] We ensure a good characterization of the intensity 
towards the limb using 285 rings, which is half size in pixels of 
the images of Fig. [4] The standard deviation of the average in- 
tensity, o"/^), was computed within each ring. As for the images 
in Fig. |U the center-to-limb variations is steeper in the case of 
the optical filter. Moreover, the fluctuations are stronger in the 
optical case with values up to 20% with respect to the average 
(Fig.S). 

We used these averaged profiles to determine LD laws for the 
simulations. For this purpose, we fitted the azimuthally averaged 
intensities determined from the simulations with a polynomial 
LD law of the form 

M=±a k{ l-pf (1) 



Table 2. Limb-darkening coefficients (see Eq.[T]with N -2) for 
the RHD simulations of Table Q] at different wavelength filters. 
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" central wavelength of the corresponding optical filter 
b central wavelength of the corresponding FLUOR filter 



as in IChiavassa et al.l (120091) . where I(p) is the intensity, a^ the 
LD coefficients, and N + I their number. We fitted all the az- 
imuthally averaged profiles using this law and weighting the fit 
by 1 /<T/(^). We varied the order N of the polynomial LD laws and 
found that N = 2 in practice already gave an optimal fit, with 
the N - 3 solution providing only minor and altogether negli- 
gible improvements to the x 2 minimization. At the same time, 
we would like to stress that using fewer terms in Eq. Q] (e.g. a 
simple power law) does not provide a good representation of the 
azimuthally averaged profiles and we therefore discourage this 
practice. In Table [2] we listed the LD coefficients a* for each 
simulation and for both filters considered in the present work. 

We have also computed classical ID, LT E, plane-parallel, 
hydrostatic MARCS mo del atmospheres dGustafsson et al.l 
1 19751: lAsplund et al.ll 1997b with identical stellar parameters, in- 
put data, and chemical compositions as the 3D simulations of 
Table Q] With Optim3D, we computed the emerging intensities 
for the same regularly spaced 285 ^-values used to get the az- 
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2.14±0.26 firn 
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x [Rsun] 

Fig. 4. Synthetic stellar disk images of the simulation with [Fe/H]=-3.0 and logg=1.6 (Table[TJ- The intensity range is 5.0xl0 3 - 
2.1 x 10 6 ergcirT 2 s _1 A -1 for the optical filter and 5.0xl0 3 -l.l x 10 5 ergcirT 2 s -1 A -1 for the FLUOR filter. 
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Fig. 5. Radially averaged intensity profiles (black line) derived from the synthetic stellar disk images of Fig. [4] and corresponding 
to RHD simulation with [Fe/H]=-3.0 and logg-\.6. The dashed grey areas denote the one sigma spatial fluctuations with respect 
to the averaged intensity profile. The intensity is normalized to the mean intensity at disk center and the radius is normalized to the 
radius given in TableQ] 



imuthally averaged intensity profiles from the synthetic stellar 
disk images and for the wavelength range spanned by both fil- 
ters, assuming a micro-turbulent broadening of 2 kms -1 . For 
each yU-value, we integrated the monochromatic intensities over 
the filters' ranges to obtain a ID intensity profile. Figure [6] il- 
lustrates the comparison between the average intensity profile 
from a 3D simulation, together with its LD fit, compared to a ID 
model profile. The latter displays large differences with respect 
to 3D profile. Moreover, full LD = fx) and partial LD 

(I(fj)/I( 1) = 0.5 ■ yu) profiles largely diverge from both ID and 
3D profiles. The effects on diameter determination from these 
visibility curves are analyzed in Section [5T2l 



5. Interferometric predictions: departure from 
spherical symmetry 



We derive the interferometric observables to determine the im- 
pact of the limb-darkening and granulation pattern on visibility 
curves and closure ph ases, using the same method described in 
IChiavassa et alj d2009h . For each synthetic stellar disk image, we 
calculated the discrete complex Fourier transform z. The visibil- 
ity, vis, is defined as the modulus \z\, of the Fourier transform 
normalized by the value of the modulus at the origin of the fre- 
quency plane, |zol, with the phase tan <p = 3(z)/!R(z), where 3(z) 
and %{z) are the imaginary and real parts of the complex num- 
ber z, respectively. We introduced also a theoretical spatial fre- 
quency scale expressed in units of inverse solar radii (Rq 1 ). The 



A. Chiavassa et al.: Three-dimensional hydrodynamical simulations of red giant stars and interferometry 

4700-5300 A 2.14±0.26 /xm 




[Fe/H] = 0.0 log(g) = 2.2 




0.6 7 



N. 

_l - ''•>- 



1.0 0.8 0.6 0.4 

4700-5300 A 



0.2 0.0 



1.0 0.8 0.6 0.4 0.2 0.0 



2.14±0.26 fim 




0.2 - 




1.0 



0.8 



0.6 



0.4 



0.2 



0.0 



1.0 



0.8 



0.6 



0.4 



0.2 



0.0 



jj, 



Fig. 6. LD fits (red dashed line) obtained using the Eq.Q]with N - 2 for the RHD azimuthally average intensity profile (solid line) of 
two simulations of TableQ] The dotted line is the intensity profile computed with the ID model having identical stellar parameters, 
input data, and chemical compositions as the 3D simulation. A full LD (dash-dotted line), a partial LD (triple dot-dashed line) are 
also shown. 



conversion between spatial frequencies expressed in the latter 
scale and in the more usual scale of 1 /" is given by: 





1 




1 


V 




= V 













d [pc] • 214.9 



(2) 



where v is the spatial frequency, 214.9 is the astronomical unit 
expressed in solar radii, and d is the distance of the observed 
star. Also useful is the following relation 



B 



v 



A ■ 0.206265 



(3) 



where v is the spatial frequency in arcsec 1 at the observed wave- 
length A in fim for the baseline B of an interferometer in meters. 

5.1. Visibility curves 

The first null point of the visibility curve is mo stly sensitive 
to the rad ial extension of the obs erved object (e.g. lOuirrenbachl 
|2001L and lChiavassa et aLll2010l for an application to RHD sim- 
ulations). In fact, the radius of Table [1] can be estimated using 



R*[Rg] 



Q[R,1 



1.22 



where vfRg 1 ] is the spatial frequency 



2-v[R-']' " L ® 

corresponding to the first null point and the angular diam- 
eter. This formula is related to the resolving power of a tele- 
scope. Our predictions concerning the first lobe of the visibility 
curves, however, cannot be used to fit the observed data because 
of the geometrical limitations of our models. In our spherical 
tiling model, based on box-in-a-star simulations, the stellar ra- 
dius is effectively a free-parameter. Only 3D simulations in the 
star-in-a-box configur ation can be strictly used for this purpose 
dChiavassa et alj2009l) . On the other hand, the first null point and 
the second lobe of the visibility curves are sensitive to the limb- 
darkening dHanburv Brown et alJll974l) . These frequency points 
are then crucial to test RHD simulations against observations to 
recover the information on the temperature stratification in the 
external layers of the star. 

Figure [7] (top panels) shows the visibility curves computed 
for 36 different cuts through the centers of the stellar disk images 
in Fig. |4] This is equivalent to generating different realizations 
of the stellar disk images with intensity maps computed for dif- 
ferent sets of randomly selected snapshots. The comparison with 



0.10 



A. Chiavassa et al.: Three-dimensional hydrodynamical simulations of red giant stars and interferometry 

4700-5300 A 2.14±0.26 fim 

0.10 



0.01 




0.01 




0.05 



0.10 



0.15 



0.20 



0.05 



0.10 



0.15 



0.20 



Ro" 



Ro" 



cm 
o 





0.05 



0.10 



Ro" 1 



0.15 



0.20 



0.05 



0.10 



0.15 



0.20 



Fig. 7. Top row: visibility curves from the RHD simulation with [Fe/H]=-3.0 and log g- 1.6 in the optical (left column) and infrared 
(right column) filters. The visibilities are computed for 36 different azimuth angles 5° apart (thin grey lines). The dotted line is an 
uniform disk scaled to match the same radius, i.e. the first null in visibility curve. The visibilities are displayed only longward of 
the first null visibility point (see text). A logarithm scale is used on y-axis. Bottom row: visibility fluctuations with respect to the 
average value as a function of spatial frequencies. 



the uniform disk visibility in top left panel of Fig. [7] shows that 
the synthetic visibilities deviate greatly from circular symmetry 
and they are systematically lower than the uniform disk model at 
the leftmost spatial frequencies shown in the plot. This is due to 
a the limb darkening effect already visible in Fig. [4] Moreover, 
bottom panels of Fig.Qshows the one <x visibility fluctuations, 
F, with respect to the average value vis (F = cr/vis): the dis- 
persion clearly increases with spatial frequency. On the top of 
the second lobe, which is easy to reach both in the optical and 
in the FLUOR filters with CHARA of a typical red giants, the 
fluctuations are about 2% (optical) and 0.3% (FLUOR). In fact, 
the effect of line blanketing in FLUOR filter (Fig.|4j right panel) 
is less important than the in the optical case and the image seems 
closer to a uniform disk. 

In addition to this, the synthetic visibilities of the simulation 
at \ogg = 1.6 and [Fe/H]=-3.0 show, in the optical filter, the 
largest deviation from UD in the range 0.17 - 0.22 R Q '. This is 
d ue to the small scale stru cture on the model stellar disk shown 
in IChiavassa et al.l d2010l) . This frequency range corresponds to 
surface structures of 2.27-2.94 Rq or 1500 -2200 Mm, approx- 
imatively the size of one of the granules as shown in Fig.Q] 



Observing red giant stars at high spatial resolution is thus 
crucial for two reasons, at least: 

(1) To further test the limb darkening predicted by our RHD 
simulations. For this purpose, it would be worthy search- 
ing for angular visibility variations, observing with the same 
telescope configuration (covering spatial frequencies higher 
than the first lobe) and using the Earth rotation in order to 
span multiple position angles. Assuming measurement er- 
rors below 2%, for visibilities of ~ 10%, one night would 
be enough. For com parison, the instrumen tal error of VEGA 
on CHARA is -1% dMourard et alj[2009b . 

(2) To characterize the granulation pattern. This requires ob- 
servations at very high spatial frequencies that may be 
possible only in the optical range. In fact, considering 
that the maximum baselin e of CHARA is 331 meters 
(Iten Brummelaar et al 1 l2005h . for a red giant with an appar - 
ent radius of 2.7 mas (e.g. HD 214868 iBaines et al.ll2010h , 
a baseline of ~270 (1200) meters is necessary to probe the 
seventh lobe at 0.5 (2. 14) jum. 

Our method of constructing realizations of red giant stellar 
disk images (Section [3} inevitably introduces some discontinu- 





Fig. 8. Top left panel: synthetic stellar disk image in the FLUOR filter for the simulation with [Fe/H] =-3.0 and log g- 1 .6 tiled using 
only one intensity map at different /i-angles (see text for details). The granulation pattern is continuous. The intensity range is [5000, 
0.11 x 10 6 ] ergcrrT 2 s _1 A -1 . Top right panel: synthetic stellar disk image in the same filter and using the same intensity map but 
flipping the y-axis every second tile to remove the periodicity at the side of each box. Bottom panel: Ratio of the fluctuations of the 
synthetic visibility curves derived from the above images: <x re f /vis le f refers to top left panel and cr shuffle /vis s h u ffie to top right panel. 



ities between neighbouring tiles by randomly selecting temporal 
snapshots and by cropping intensity maps at high latitudes and 
longitudes. We checked whether these artifacts produce signif- 
icant fluctuations in the theoretical visibility curves. These dis- 
continuities in the granulation pattern stem from using random 
snapshots of the simulation sequences. We explored the impact 
of such discontinuities at the edges of neighbouring tiles by con- 
structing two mock stellar disk images (Fig. [8j. To generate the 
first image, we tiled a sphere with intensity maps computed for 
one single simulation snapshot: in this case, the continuity of 
the granulation pattern is ensured by the the periodic horizontal 
boundary conditions of the simulation. For the second image, we 
used the same intensity map as for the previous case but flipping 
the y-axis every second tile: in this way, we generated a gran- 
ulation pattern that is statistically similar to the one in the first 
image in terms of fine spatial structures but at the same time has 
discontinuities at the boundaries between tiles (top right panel). 
We then computed the visibility curves for 36 different angles 
5° apart and derived the visibility fluctuations (Fig. [8] bottom 
panel). Both in the optical and in the FLUOR filter, the visibility 
fluctuations caused by the disentangled periodicity at the sides 



of each boxes have a contribution about two orders of magnitude 
smaller than the signal observed in Fig. [7] We conclude that our 
approach to construct the orthographic projection with random 
snapshots have a smaller effect on the visibility curves than the 
signal originated by the inhomogeneities of the stellar surface. 

In relatively large filters, such as FLUOR, several spatial 
frequencies are simultaneously observed by th e interferome- 
ter. Th is effect is called bandwidth smearing. iKervella et al.l 
show that this effect is negligible for squared visibili- 
ties larger than 40% but it is important spatial frequencies close 
to the first minimum of the visibil ity function. To account for 
this effect, IWittkowski et al.l ( |2004|) proposed to compute aver- 
aged squared visibility amplitudes, (vis 2 ), as 



(vis 2 ) 



I' vis l dA 



JG 



vll 



TjFjdA 



(4) 



where vis 2 is the squared visibility at wavelength A, Aq and 
A\ are the filter wavelength limits, Tx is the transmission curve 
of the filter and F\ is the flux at wavelength A (see Fig. [3). 
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Fig. 9. Visibility fluctuations with respect to the average value 
as a function of spatial frequencies for the RHD simulation with 
[Fe/H]=0.0 and logg=2.2 in the FLUOR filter. The dashed line 
corresponds to fluctuations of visibility computed taking into ac- 
count the bandwidth smearing effect (see text and Eq. |4]i while 
the solid line to fluctuations of visibility computed without band- 
width smearing. 



In our case, first we generated a synthetic stellar disk image at 
every wavelength and then we computed the visibility curves 
vis\ for 36 different cuts through the centers of the stellar disk 
images. Then we applied Eq. [4] to obtain the averaged squared 
visibilities and finally we derived the visibility fluctuations with 
respect to the average values. Figure [9] shows that there is not 
a clear difference between the visibility fluctuations computed 
with and without taking into account the bandwidth smearing 
effect. 

The simulations of ICollet et all (l2007l) showed that the pre- 
dicted characteristic size of granules decreases with the metal- 
licity of the star. Therefore, we focus now on the simulations of 
Table [TJ with the same surface gravity, log g-2.2, and different 
[Fe/H] and effective temperature. We considered only the optical 
filter that show larger (and thus easier to detect) visibility fluctua- 
tions. Figure[10]displays the synthetic stellar disk images and the 
visibility fluctuations derived from the visibility curves. The so- 
lar metallicity simulation has almost systematically smaller fluc- 
tuations (except for the 3rd lobe at R ' ~ 0.10) at all frequencies 
until Rq 1 ~ 0.26 (i.e. the 6th lobe, that is close to the baseline ex- 
tension limit of CHARA interferometer for a 2.7 mas star). The 
difference between [Fe/H]=0.0 (black curve) and [Fe/H]=-3.0 
(light blue curve) is about 1 % (same order of the instrumental 
error of VEGA) on the top of the second lobe (R ' ~ 0.06) with 
peaks of about 1.2-1.7 % at R Q ' ~ 0.17 and 0.20, respectively. 
We emphasize that it is possible to study the metallicity depen- 
dence on visibility curves to characterize the granulation pat- 
tern of red giant stars. It is however important to carefully con- 
sider specific instrumental effects that could affect the visibility 
fluctuations. This can be already done with p resent day interfer- 
ometers in addition to image reconstruction dBerger et al.ll2010l 
iMalbet et ai1l2010l) which eventually will be able to measure the 
size of granules directly. 



5.2. Estimation of diameter correction: 1D versus 3D models 

Interferometric data are often fitted with a simple uniform disk 
model that, although unphysical, has the advantage of immedi- 
ately telling if a star is resolved, after which limb-darkening cor- 
rections are applied. Thus, to estimate the angular diameter from 
measured visibilities it is necessary to know the intensity distri- 
bution of the stellar disk and, for this purpose, parametric (e.g. 
full or partial limb darkened) or photospheric one dimensional 
models (e.g. MARCS) models are used. 

We used the intensity profiles obtained from 3D and ID 
models (see Section |4] and Fig.|6]l to estimate the correction on 
the uniform disk diameter determination. We derived the visi- 
bility law Va (B,®) from the intensity profile I(A,fi) using the 
Hankel integral: 



V A (B,®) = 




where A is the wavelength in meters (in this case the central 
wavelength of the filters in Fig. |3), B is the baseline in meters, 
is an arbitrary angular diameter in radians (we assume 2 mas), 
/o the zeroth order of the Bessel function, fi - cos(6) as defined 
in Section|4] and A the normalization factor: 

A = f I(A,[i)ndn (6) 
Jo 

Figure [TTJ shows the visibilities computed from the inten- 
sity profiles of Fig. [6] For typical red giant stars having 4600 < 
T e ff < 5100 K (compare with Table [1} differences in angular 
diameters vary from ~ -3.5% to ~ 1% in the optical, and are 
roughly in the range between -0.5 and -1.5% in the infrared 
(TableO, the corresponding change in effective temperature be- 
ing Ar e ff/r e ff = 1 - V®3d/®id- These corrections depend on 
wavelength, the deviations being moderately stronger in the op- 
tical than in the infrared. Differences between the ID and the 
mean 3D atmospheric structures and inhomogeneities of red gi- 
ant stars have been already studied by means of 3D hydrody- 
namical simulations dCollet et al .1 12007). These differences have 
a significant effects on the predicted strengths of spectral lines as 
well as on the center-to-limb profiles. As a con sequence, wh en 
using LD laws obtained from ID models (e.g. IClarell2000l) to 
determine angular diameters from observed data, one must take 
into account that 3D models usually predict a lower center-to- 
limb variations with respect to ID, thus implying smaller radii 
and moderately hotter effective temperatures. 

While the impact of the corrections in Table[3]is usually not 
dramatic, they are not negligible to properly set the zero point 
of the effective temperature scale derived by mean of this fun- 
damental method. In particular, when future long baseline inter- 
ferometers working in the optical will be able to (partly) resolve 
(very) metal-poor stars, it will be fundamental to take those cor- 
rections into account for a correct derivation of their diameters. 
Moreover, the model dependence on the diameter determination 
has also an impact on the reliability of the existing catalogues 
of calibrator stars for interferometry that are based on red giants 
with similar parameters to those studied in this work. In fact, 
the formal hig h accuracy reached on diameter determination (a 
fraction of 1%, Me rand et al.ll2010T) is based only on ID models, 
while we show that the correction reported in Table 0c an have a 
sizeable impact. 

Our results i n Table |3] are in qualitative agreement with 
those o btai ned bvlAllende Prieto et alj d2002l) . lAufdenberg et al.l 
(l2005h and lBigot et al.1 d2006l) for dwarf stars, though the lower 
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Fig. 10. Top and central rows: synthetic stellar disk images in the optical filter for the four simulations of Table [1] having log g-2.2. 
The intensity range is [5000, 1.5 x 10 6 ] erg cirT 2 s" 1 A" 1 for [Fe/H] = 0.0 and -1.0; [5000, 2.0 x 10 6 ] erg cirT 2 s" 1 A" 1 for [Fe/H] = 
-2.0 and -3.0. Bottom panel: fluctuations of the synthetic visibility curves derived from the above images. The colors have the same 
meaning as in Fig. [3] We dropped the first lobe of the visibility curves as in Fig. [7] 



surface gravity of red giants might also play a role, since the 
size of surface related convective structures is roughly inversely 
proportional to the gravity. As we already discussed, the size of 
these inhomogeneities have an impact on the average intensity 
profile and thus on the interferometric observables. 



These corrections could also be relevant in explaining some 
discrepancies known in literature, which involve the compar- 
ison with interferometric angular diamet ers. For example, i n 
the library of stellar spectra assembled by ICohen et al.l (Il999l) . 
absolute spectrophotometry has the tendency to return angu- 
lar diameters systematically smaller by a few percent with re- 
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Fig. 11. Visibility curves computed with Eq.[5]for and angular diameter of 2 mas and for the 3D (solid line), ID (dotted line), and 
fully limb darkening (dash-dotted line) intensity profiles of Fig. [6] In top panels, the dotted and solid lines are almost overlying, 
while in bottom left panel the dotted line is almost overlying the dash-dotted one. A logarithm scale is used on y-axis. Synthetic 
visibilities in these plots are not realistic near the nulls but are intended only for model-to-model comparison. 



Table 3. Ratio between the limb-darkened diameters recovered 
using ID MARCS models (©id) or 3D simulations (©3d) and 
the corresponding change in effective temperature Ar e jf for the 
RHD simulations of Table Q] 



A \pm] 


[Fe/H] 


logg 


®3d/©1D 


Ar eff [K] 


0.5^ 


0.0 


2.2 


1.003 


-7 


2.14~ 


0.0 


2.2 


0.996 


9 


0.5 


-1.0 


2.2 


0.991 


21 


2.14 


-1.0 


2.2 


0.996 


9 


0.5 


-2.0 


2.2 


0.982 


46 


2.14 


-2.0 


2.2 


0.990 


25 


0.5 


-3.0 


2.2 


1.011 


-28 


2.14 


-3.0 


2.2 


0.989 


28 


0.5 


-3.0 


1.6 


0.965 


82 


2.14 


-3.0 


1.6 


0.984 


37 



" central wavelength of the corresponding optical filter 
* central wavelength of the corresponding FLUOR filter 



spect to interferometry. Smaller angular diameter s for cool gi- 
ants w ould also explain the dichotomy reported bv lAlonso et al.l 
(1 19941) when trying to set on the same scale r e ff derived for hot 



and cool stars from interferometry with those determined from 
the Infrared Flux Method. Recent revision of the latter technique 
support in fact hotter effective temperatu res also for dwarfs and 
subgiants dCasagrande et al.l2006ll20Toh and the extension to gi- 
ants is currently under way (Casagrande et al. in prep.). Smaller 
limb-darkened angular diameters (i.e. hotter effective temper- 
atures) would also reduce the discrepancy recently noticed by 
Bain es et al 1 (120101) when comparing interferometric with spec- 
troscopic T e ff of giants. Thus, it is important to keep in mind 
that even when interferometry is used to derive effective temper- 
atures in the most fundamental way, there is a non negligible - 
and still partly unexplored- model dependence in the correction 
from uniform to limb-darkened disks. We remark that our results 
are currently limited to a rather small range in effective temper- 
atures and gravities and to this respect, 3D models will also give 
an important contribution to set the effective temperature scale 
on a physically sounder base. 

5.3. Closure phases 

One can use closure phase between three telescopes, as the sum 
of all phase differences removes the atmospheric contribution, 
leaving the phase information of the object visibility unaltered 
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Fig. 12. Scatter plot of closure phases of 2000 random baseline 
triangles for RHD simulation with [Fe/H]=-3.0 and logg=1.6 
with a maximum linear extension of ~510 m at 0.5 fim and 
~2200 m at 2.14 /mi and for a star at 100 pc. The extension 
of the baselines has been expressed in R Q ' using Eq. [2] and [3] 
The axisymmetric case is represented by the straigth lines. 



(e.g. lMonn ier 2007). The closure phase thus offers an important 
complementary piece of information, revealing asymmetries and 
inhomogeneities of stellar disk images. 

Figure [T2l displays large deviations from axisymmetric case 
(zero or +n). This is an evident signature of surface inhomogen- 
ities. There is a correlation between Fig. [12] and Fig.|7]because 
the scatter of closure phases increase with spatial frequencies as 
for the visibilities. Moreover, the optical filter shows larger de- 
viations than the FLUOR one, and in both cases there are asym- 
metries already detectable at frequencies corresponding to the 
top of the second lobe. It would be very efficient to constrain 
the level of asymmetry and inhomogeneity of stellar disk im- 
ages by accumulating statistics on closure phase at short and 
long baselines. Small deviations from zero will immediately re- 
veal departure from symmetry thus complementing information 
of the visibility curves towards a better characterization of limb 
darkening and granulation patterns. The CHARA interferometer 
is well adapted for this purpose because it can combine several 
telescope at once and get reliable closure phases. 

The characterization of the closure phase, together with 
the limb darkening behavior (see Section [Hi, is important for 
two more reasons in red giant stars. First, iRenard et a il (120081) 



demonstrated that closure phases provide the differential planet 
to star contrast ratio as a function of wavelength. When observ- 
ing a star with a faint companion, their fringe patterns add up in- 
coherently and the presence of a planet causes a slight decrease 
in the phase changes and, consequently, the closure phases. This 
difference can be measured with a temporal survey and should 
be corre cted with the intr insic closure phases of the parental star. 
Second. rvan Belld (|2008) provided a technique based on closure 
phases to determine the orbital plane position angle of a planet 
transiting in front of the star as well as the planet's radius us- 
ing the high-precisio n multi-telescope beam combiner (MIRC) 
dMonnier et al.ll2006b on CHARA. Also in this case, theoretical 
predictions of the closure phases of the parental stars are crucial. 

6. Comparison with the interferometric 
observations of HD 214868 

In this Section, we compare our RHD simulations to the gi- 
ant HD 214868 (K2 III, r eff = 4440 + 70 logg = 2.1 + 0.2 
and [Fe/H]= -0.18 + 0.50 dex: lBaines et~aT1l2O10l) . This star has 
been chosen because it has been observed up to the second lobe 
with CHARA in a filter cent ered at ~ 2.15 nm an d sensitive 
to the range 1.00 - 2.50 /im dMcAlister et alJl2005h . which is 
close to FLUOR filter. For this purpose, we chose the model of 
Table Q] having the closest parameters (i.e. T e ff = 4697 + 18, 
[Fe/H]=0.0, and log g=2.2) and compared our synthetic visi- 
bilities in the FLUOR filter to the observations (Fig. [T3l . The 
absolute model dimensions have been scaled to match the in- 
terferometric observation in the first lobe. However, as already 
mentioned in Section [5] the information about the radius can- 
not be retrieved because of the limitations in the method we 
use to construct the stellar disk images. We therefore concen- 
trate on the visibility data points in the s econd lobe (i.e. where 
the li mb darkening effects are important; lHanburv Brown et al.l 
Il974l) . A Levenberg-Marquardt least-squares minimization re- 
turns reduced ;f 2 =0. 59 for the best RHD synthetic visibility, all 
RHD visibilities falling in the range [0.59, 0.66]. Formally, the 
RHD simulation shows an improvement over parametric model 
(uniform disk model has reduced^- 2 =0.72) for the interpretation 
of these interferometric observations, though it must be noted 
that observed points in Fig. [13] fall at baselines where 3D and 
uniform disk models nearly overlap. A more uniform sampling 
(e.g. close to the first null point) will be very valuable to fur- 
ther test our simulation. For the sake of completeness, we also 
remark that the simulation used in this specific exercise is not 
tailored exactly to the stellar parameters of the studied star, and 
improvements along this direction are still possible. 

7. Conclusions 

In this work, we provided a set of interferometric predictions 
from our RHD simulations of red giant stars that will be tested 
against observations with today interferometers like CHARA. 
Red giant stars are ideal targets for interferometers because they 
are numerous and have large angular diameters as well as bright 
infrared apparent magnitudes. 

We derived average center-to-limb intensity profiles from 
the synthetic stellar disk images based on RHD simulations and 
computed limb-darkening coefficients in the optical as well as in 
the infrared. The predicted center-to-limb variations of 3D sim- 
ulations can be tested against observations: the visibilities and 
the closure phases show evident departures from from circular 
symmetry, which are due to inhomogeneities on the stellar sur- 
face. We emphasize that red giant stars should be observed at 
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high spatial frequencies, to characterize the granulation pattern, 
and in particular in the optical range where the visibility fluctua- 
tions are larger. In addition to this, it would be valuable to couple 
the observations in the optical and in the infrared to test models 
predictions for limb darkening. We counted at l east 58 red gi- 
ants st ars using the stellar diameters catalogue of lLafrasse et al.l 
(120101) with stellar parameters corresponding to our simulations, 
an apparent diameter larger than 2 mas in the visible and V mag- 
nitude lower than 6 (the current magnitude limit of VEGA is 7, 
iMourard et al.|[2Q09h . We provided also a first comparison with 
the red giant star HD 214868 showing that our RHD simula- 
tions returns a better fit with respect to simple parametric uni- 
form disk models. More observations along this direction will 
provide useful information not only on red giants, but also on 
our RHD models. This will include the center-to-limb variation, 
which carry the information on the temperature stratification in 
the stellar atmosphere, and the characterization of the granula- 
tion pattern in terms of size and contrast with varying metallic- 
ity and effective temperature. Moreover, the study of the closure 
phases behavior will help to recover the physical parameters of 
possible hot Jupiters around red giant stars. 

Interferometry is advantageous in that it provides the abil- 
ity to directly measure stellar angular diameters. We found that, 
for typical red giant stars studied in this work, differences in an- 
gular diameters using 3D or ID limb-darkening laws vary from 
~ -3.5% to ~ 1% in the optical, and are roughly in the range 
between -0.5 and -1 .5% in the infrared, with 3D models return- 
ing in general smaller angular diameters with respect to ID, thus 
implying moderately hotter effective temperatures. The impact 
of these corrections is not dramatic, but they are not negligible 
to correctly set the zero point of the effective temperature scale 
derived by means of this fundamental method, in particular, for 
future observations of metal-poor stars. Moreover, the model de- 
pendence on the diameter determination has also an impact on 
the reliability of the existing catalogues of calibrator stars for 
interferometry that are based on giants with parameters similar 
to those covered in this work. We emphasize that even when in- 
terferometry is used to derive effective temperatures in the most 
fundamental way, there is a non-negligible and still partly unex- 
plored model dependence. 
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Fig. 13. Comparison of the RHD simulation with [Fe/H]=0. and 
log g-2.2 (Table [T|) to the observa tions (red dots) of the gi- 
ant HD 214868 dBaines et al.ll2010h . The solid thick line cor- 
responds to the best fit (reduced ;^ 2 =0.59) visibility curve, 
the dashed line is the unifor m disk (2.72 1 mas with reduced 
X 2 =0.72) model introduced bv lBaines et al.l . The visibility curve 
from RHD simulation never goes to zero due to the bandwidth 
smearing effect. A logarithm scale is used on y-axis. 
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